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ABSTRACT 

Context. The effective temperature scale of FGK stars, especially at the lowest metallicities remains a major problem in the chemical 
abundance analysis of metal-poor stars. 

Aims. We present a new implementation of the infrared flux method (IRFM) using the 2MASS catalogue. 
' Methods. We computed the theoretical quantities in the 2MASS JHK S filters by integrating theoretical fluxes computed from ATLAS 

models, and compare them directly with the observed 2MASS JHK S magnitudes. This is the main difference between our imple- 
mentation of the IRFM and that of Ramirez & Melendez (2005, ApJ, 626, 446; hereafter RM05), since to introduce new stars at the 
lowest metallicities they transform the 2MASS JHK^ magnitudes into the TCS photometric system. We merge in our sample the 
stars from Alonso et al. (1996, A&AS, 117, 227; hereafter AAM96; 1999, A&AS, 139, 335; hereafter AAM99), and other studies 
to appropriately cover a wide range of metallicities, ending up with 555 dwarf and subgiant field stars and 264 giant field stars. We 
derived a new bolometric flux calibration using the available Johnson-Cousins UBV(RI)c and the 2MASS JHK S photometry. We also 
computed new r cff versus colour empirical calibrations using our extended sample of stars. 
Results. We derived effectives temperatures for almost all the stars in the AAM96 and AAM99 samples and find that our scales of 
temperature are hotter by ~ 64 K (cr T = 104 K, N = 332 dwarfs) and ~ 54 K with a <r T = 131 K (N = 202 giants). The same 

comparison with the sample of RM05 for stars with [Fe/H] < -2.5 provides a difference of 87 K (cr T = 194 K, N = 12 dwarf 

, stars) and ~ 61 K (cr r = 62 K, N = 18 giant stars). 

Conclusions. Our temperature scale is slightly hotter than that of AAM96 and RM05 for metal-rich dwarf stars but cooler than that 
' of RM05 for metal-poor dwarfs. We have performed an fully self-consistent IRFM in the 2MASS photometric system. For those 

' who wish to use 2MASS photometry and colour-temperature calibrations to derive effective temperatures, especially for metal-poor 

stars, we recommend our calibrations over others available in the literature. In our implementation we avoid the transformation of the 
. 2MASS JHK S magnitudes to a different photometric system and thus fully exploit the excellent internal consistency of the 2MASS 

photometric system. 

Key words, infrared: stars - stars: abundances - stars: atmospheres - stars: fundamental parameters 

1 . Introduction than on a direct measure of th e angular diameter. Th e infra red 

flux method (hereafter IRFM; Black welTet al.| [l990. and refer- 
The effective temperature is a function of the bolometric flux and encefj therein) ifj eciall ad te for determining the e ffec- 

g ; the angular diameter according to the equation tive temperature o f F , G and K stars. The IRFM was first in- 

4 troduced bv lBlackwell & Shallisl (I1977I) who proposed simulta- 

T e (f = (—)^ 4 6 1/,2 ^boi (1) neously determining the effective temperature and the angular 

°~ diameter of a star. The basic idea is to use the monochromatic 

where cr is the Stefan-Boltzmann constant, 9 the angular diame- fl ux m the infrared since it is mainly dependent on the angular 

ter, and F bo i the bolometric flux measured on the surface of the diameter but is approximately dependent only on the first power 

Earth. However, direct measure of angular diamet ers is restricted f r eS , whereas the integrated flux strongly depends on the tent- 



er 
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to relatively few stars, especially for dwarf sta rs. Kervella et alj perature (proportional to T 4 ff ) 



(2004, [20081 interferometry) and lBrown et alj d200U transit ob 

servations) have directly measured the angula r diameters of Popular indirect methods for deriving effective temperatures 



bright stars. Only recently, Bai n^tifl <|2008) have used the areflieexcitation equilibrium of Fei lines (e.g. |Santos et al 



CHARA interferometric array to provide measurements of an- 2004 ' 200 ^ and on fitting Balmer lines (e.g. | Fuhrmannet al 



gular diameters of -28 dwarf and subgiant stars, although all of ll994j; | Barklem et al. | |2002j). Temperatures based on 

them have metallicities [Fe/H] > -0 5 ^ e 1 exc i tat i° n equilibrium depend on the model assumptions, 

A semi-direct method of temperature determination is one such as non-LTE effects, especially in metal-poor stersjsee 

that makes use of Eq. H]but relies on model atmospheres, rather iTfievenin & Idiart| |1999t | Shchukina & Truiillo Bueno | |2001j). 

Recently, Barklem (2007) has also raised concerns about pos- 

Send offprint requests to: J. I. Gonzalez Hernandez sible non-LTE effects on the w ings of Balmer line s. Both excita- 

* Cosmological Impact of the First STars (CIFIST) Marie Curie tion e quilibria dAsplundl l 2005l) and Balmer lines dLudwig et all 

Excellence Team I2009L in prep.) are also sensitive to granulation effects. This 
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makes such methods strongly model-dependent, which is an un- 
desirable feature. However, temperatures derived from Balmer 
lines and Fe i excitation equilibria have the considerable advan- 
tage of being reddening independent. 

One of the motivations of this work is to investigate the trend 
of Li abundances towards low metallicities ([ Fe/H] < —2.5), 
using our own implementation of the IRFM. dBonifacio et alj 
2007) investigated the Spite plateau at the lowest metallici- 
ties (down to [Fe/H]=-3.3) and found marginal evidence that 
there could be an increased scatter or even a sharp drop in the 
Li abundance at these low metallicities. Determination of the 
baryonic density from the fluctuations in th e cosmic microwave 
backg round (CMB) by the WMAP satellite dSpergel et al.ll2003l 
|2007|) implies a primordial Li abundance, which is at least 
a factor of 3-4 larger than observed on the Spite plateau, 
creatin g a conflict w i th the traditional interpretation of the 
plateau dSpite & Spite! Il982allbl) . This discrepancy would be 
even greater if the drop i n the Li abundance v ersus metallicity 
were to be confirmed (see Sbordone et al. 2008, in prep.). 



One decade ago, iBonifacio & M olar 3 (119971) determined 
Li ab undances using the IRFM temperatures of lAlonso et al.l 
d!996al) . They investigated the different Li abundance trends 
found with different temperature scales. In particular, the tem- 

which are based on the IRFM 
are cooler than the IRFM 



nperature 
(fl996h . 



peratures of iRvan et al.l 
implementation of Magain (1987), 
temperatures of lAlonso et al.l(ll996al) . provided r A ionso - ^Ryan ~ 
+ 10 K at [Fe/H] ~ -1.7 and ~ + 130 K at [Fe/H] ~ -3.3. From 
this, [Bonifacio & Molarol d 19971) concluded that the presence or 
absence of trends in lithium abundance with T e ff is strongly de- 
pendent on the temperature scale adopted. 

Mele ndez & Ra mirez (2004) applied their own IRFM imple- 
mentation ( Ramirez & Melendez 2005a) to deriving the effec- 
tive tempe rature and Li abun dances for a sample of stars similar 
to that of Rya n et all {1996). They find individual temperature 
differences of up to 400-500 K for the some stars with metal- 
licity below -3.0 dex. More recently, IBonifacio et al.l ([2007) 
have compared the temperatures obtained from Ha profiles to 
other temperature indicators, among them those from the IRFM- 
based colour-t emperature c alibrati ons of dRamfrez & Melendez! 
2005a) and lAlonso et al. Ul996b). When a reddening based 
on the ISchlegel et al.l d 19981) maps is ado pted, from the 
r eff :(V-K) calibration of lAlonso etail d!996bl) . the mean dif- 
ference 7V-a:) aam% - T Ho is only 8K with a standard devia- 
tion of 100 K. However, if we use the T e ff:(V - K) calibra- 
tion of Ramirez & Melendez! d200 5b). this mean difference is 
7\v-a")rmo5 _ T Ha of 265 K, with a standard deviation of 122K. 

Ramf rez & Melendez! (|2005a) add as calibrators a small 
samp le of metal-poor sta rs mainly from IChristlieb et ail d2004l) 
and ICavrel et al. (2004), and a larger sample of metal-rich 
stars from ISantos et al.l d2004l) to the original sample of 
lAlonso et al.1 d!996al) . and computed new r e jf-colour calibra- 
tions. Since the majority of calibrators shared by the two sam- 
pl es, this large differenc e (~ 25 K) between the c alibr ations 
of iRamfrez & Melendez! d2005al) and lAlonso et all d!996bl) at 
low metallicity is a priori unexpected. One could argue that 
the models used by the two groups are not exactly the same, 
however they must be very similar (ATLAS 9 models with the 
same ODFs and microturbulent velocity). Since the IRFM is 
only wea kly dependent on the mod els adopted, as shown by the 
results of ICasagrande et all (|2006), who used both ATLAS and 
MARCS models, it seems unlikely that this difference is rooted 
in the different models. We suggest instead that this is because 




A(nm) 



Fig. 1. Spectral energy distribution of the ATLAS model of Vega 
(T eS = 9550 K, logg = 3.95, [Fe/H] = -0.5 and v micl , = 
2 km s 1 ) in the infrared. The solid and dashed-dotted lines rep- 
resent the intensity and continuum flux, respective ly. The filled 
circles are the monochromatic fluxes adopted by ICohen et al.l 
d2003l) . We also show the transmission functions of the JHK S 
2MASS filters. 



IRamrrez & Melendez! {2005a) use 2MASSQ JHK S magnitudes 
for the low metallicity calibrators; such magnitudes were then 
transformed into the TCS system to merge them with th e ho- 
mogeneous set of TCS photometry of Alons o et al.l (fl996a). The 
errors in the transformation between the 2MASS and the TCS 
systems are then added to the photometric error and may have 
undesired effects on the final calibration. We have no way to 
prove that this is indeed the case; however, to circumvent such 
problems in this paper, we propo se a new implementation of 
the IRFM, including the stars fromlAlonso et al.1 d!996al Il999al) 
and IRamrrez & Melendez (2005a), but using the 2MASS pho- 
tometry for all calibrators rather than a mixture of 2MASS and 
TCS. The 2MASS magnitudes are probably not as accurate as 
the careful TCS photometry of lAlonso et alj dl996al) . but the in- 
ternal consistency of the 2MASS photometry is very high, about 
1-2% dCutri et alJl2003l) . 



2. Implementation of the IRFM 

The IRFM dBlackwell et alj 1 19901) evaluates the quotient be- 
tween the bolometric flux, Fboi, and the monochromatic flux at 
a chosen infrared wavelength, F (Am), both measured at the sur- 
face of the earth, as an indicator of the T e g . This quotient is the 
so-called observational 7?-factor, R a \, s . The theoretical counter- 
part derived from models, R t h eo , is obtained as the quotient be- 
tween the integrated flux, 0"!Tf ff , and the monochromatic flux at 
^ir, ^mod(^iR), at the surface of the star. Thus the basic equation 
of the IRFM is 



^obs - 



Fbo 



F (/Iir) F mod (Ai R ,T eS , [Fe/H], log g) 
= Rtheo(A]R,T eS , [Fe/H], log g) 



(2) 



where the dependence of models on metallicity, surface grav- 
ity, and /Iir is explicitly taken into account. The monochromatic 



1 The 2MASS c atalogue can 
http://www.ipac.caltech.edu/2mass/ 



be 



accessed at 
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fluxes are obtained by applying the relation 

F(Am) = q(A m ,T eS , [Fe/H], logg)[ J F cal (^ IR )10-°- 4(m *-" , -' ) ] (3) 

where ra+ is the magnitude of the target star, and m cd \ and F c - d i 
are, respectively, the magnitudes and the absolute monochro- 
matic fluxes of the calibrator star (see Table[T]and Sect. [5]). The 
^-factor, usually ~ 1, is a dimensionless factor that corrects the 
effect of the differe nt curvature of the flux distribution across the 
filter window (see lAlonso et alj|1994 Il996al Il999al for more 
details). We have used the definition of Alo nso et al.l d!996al) for 
the computations of the ^-factors (see Sect. [5]). 

By merging the previous two equations we can separate the 
observational and model inputs as 



fbol 



^calWlR)10- - 4 ('"*- m ">> 

qiMJeff, [Fe/H], logg)R theo (AiRJ e n; [Fe/H], logg) (4) 

The synthetic magnitudes, the q- and R- factors, necessary for 
implementing of the IRFM were computed f rom the ATLAS the- 
oretical fluxes of Castelli & Kurucz (20030 using the 2MASS 
JHK S filters. We used the fluxes in the ranges 3500 K< T e ff < 
7500 K, 0.0 <logg < 5.0, and -4 < [Fe/H] < +0.5, and for 
metal-poor models with [Fe/H] < -0.5, we used the fluxes com- 
puted from the a-enhanced models. We derived a new calibra- 
tion of the bolometric flux in the 2MASS photometric system 
(see Sect. |6}. 

3. Online data available at the CDS 

Several tables are available at the CDS0. We provide eight tables 
containing the q- and /^-factors computed as indicated in Sect. [2] 
and for metallicities in the range [^4,0.5], temperatures in 
the range [3500,50000], and gravities in the range [0,5]. Within 
these tables, we also put the theoretical colour V-K and magni- 
tudes JHKs in the 2MASS system for each atmospheric model. 
These theoretical colour and magnitudes, which are not used in 
this work, weren normalised to Vega assuming V=J=H=Ks=0. 
If the user wants to use a different zero point for Vega, it is trivial 
to add it to our theoretical magnitudes. In addition, eight tables 
containing the A -factors and BX-B coefficients for the same set 
of models are also available at the CDS, needed for the bolomet- 
ric flux calibration (see Sect.|6]l. At the CDS, we also provide two 
tables, with 555 dwarf stars and with 264 giant stars, containing 
the photometric data and reddenings used in this paper, stellar 
parameters and metallicity, bolometric fluxes and IRFM temper- 
atures for all the stars in our samples (see Sect. 14. Il and l4.2b . 

4. Sample, observational data, and stellar 
parameters 

4.1. Photometric data 



Our sa mple includes almost all the stars in lAlonso et al. (1996a, 
1999a) with available photometric data in th e final release of 
the 2MASS catalogue (ISkrutskie et al.l 120061) and with photo- 
metric accuracy < 0.3 mag, for temperature determinations. 
We adopted th i s rathe r high tolerance because giant stars 
of lAlonso et al.l (fl99 9a) are relatively bright and usually the 
2MASS photometric accuracy is very low for these stars. 



However, to improve the precision of the bolometric flux and 
r e ff:colour-[Fe/H] calibrations significantly, we decided to fur- 
ther lower the accuracy limit down to < 0.1 mag (see Sect. [6] 
and [9]). Therefore, stars with 2MASS photometric errors > 0.1 
mag were only used for the purpose of deriving effective tem- 
peratures and they are provided as online data at the CDS. 

We adopted the same UBV(RI) c photometric data as used 
by iRamfrez & Melendezl d2005albl) . which were kindly pro- 
vided by Ramirez & Melendez (private communication). These 
data were extracted from the G eneral Catalogue of Photometric 
Data dMerrnilliodet a"j]|l997l GCPD). For those stars of the 
lAlonso et al.l d!996al Il999al) samples without V data in the 
GCPD these magnitudes were obtained from SimbacQ, and were 
later used to derive bolometric fluxes and IRFM temperatures. 

4.2. Metallicity and surface gravity 



For our sample of stars from (lAlonso et al. | Il996al [[999a), 
we adopted the surface gravities and metallicities provided by 
Ramirez & Melendez (private commu nication) which mostly 
use th e mean values of those reported in lCavrel de Strobel et al.l 
d200lh . 

We completed our sample of dwarfs and subgiants with 
the metal- rich and metal-poor s tars a lready included in the 
sample of IRamfrez & Melendezl d2005al) . The metal-rich sam- 
ple mostly contains planet-host stars and the comparison 
sample from ISantos et"aT] d2004l) . but we a lso added to our 
sampl e the stars with [Fe/H] > -2 from ICasagrande et al.l 
(2006). We completed the s ample with extremely metal -poor 
dwarfs fro m Bonifa cio etail (120071) . IChristlieb etall (|2004|) . and 
Barklem et al. (2005). For these stars, we adopted the same sur- 
face gravity and metallicity as published in the abo ve papers. 



2 http://wwwuser.oats.inaf.it/castelli/grids.html 

3 http://cdsweb.u-strasbg.fr/ 



Our sample of giants contains the stars in lA lonso et al 
(1999a), plus the metal-poor stars from IRamfrez & Melendez 
(2005a ). This includes stars from the "First Stars" project 
(ICavrel et all 12004 ISoite et al.ll2005h . and we adopted the sur- 
face gravity and metallicity for these stars as provided in these 
papers. 

The errors on surface gravity and metallicity for all dwarf, 
subgiant, and giant stars were assumed to be A logg = 0.5 dex 
and A[Fe/H] = 0.1. The average systematic errors due to a differ- 
ent metallicity (by +0.1 dex) and a different surface gravity (by 
+0.5 dex) are 13 K and 1 1 K, respectively, for dwarfs, and 1 1 K 
and 28 K for giants. These errors were estimated by quadrati- 
cally adding the errors on effective temperature from each band 
and calculating the average over all stars in both samples. 

4.3. Reddening corrections 

The extinction in each photometric band, A,-, as determined us- 
ing the relation A, = RjE(B — V), w here is given by the 
coefficients provided in (McCall 2004). Reddening corrections, 
E(B-V), were adopted from Ramirez & Melendez (private com- 
munication). For the metal-rich stars of Casagra nde et al.1 (12006) 
and the extremely metal- poor dwarfs of iBonifacio et al.l (120071) 
and IChristlieb et all (2004), redd ening corrections were derived 
from the maps of ISchlegel eTail (1 19981) . The E(B - V) from the 
maps is appropriate for objects outside the dust layer, which 
is confined to the Galactic disc. For objects which are within 
the dust layer the map estimate should be corrected by a fac- 
tor [1 - exp(-|<isinfo|//i)], where d is the distance of the star, 
b its galactic latitude and h the scale height of the dust layer 



http://simbad.u-strasbg.fr/simbad/ 
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Fig. 2. Left: Empirical calibration cf>:(V - K s )-{Fe/H] for dwarfs in the metallicity bins -0.5 < [Fe/H] < 0.5 (filled circles), -1.5 < 
[Fe/H] < -0.5 {open circles), -2.5 < [Fe/H] < -1.5 (triangles), and [Fe/H] < -2.5 (diamonds). The lines correspond to our 
calibration for [Fe/H] = (solid line), -1.0 (dotted line), -2.0 (dashed line), -3.0 (dotted-dashed line). Right: Residuals of the fit 
(Ac/> = (cp c . d i - <^irfm)Mrfm) as a function of (V - K s ) and [Fe/H]. 



(see, e.g. [Bonifacio et al. 2000b). For this purpose we used the 
parallaxes provided by Simbad (which co me mainly from the 
Hipparcos catalogue lPerrvman et al.ll 19971) and assumed a scale 
height of the dust laye r of 125 pc. Bonifaci o et alj (HbOOa) note 
that, when the maps of Schle gel et alJ dT998) provide reddenings 
larger than 0.1 mag, they overestimate the reddening with re- 
spect to other indicators, and proposed a simple formula for cor- 
rectin g the reddening from th e maps. We make use of formula 
(1) of Bonifaci oetaLl d2 000a: i to corre ct the reddenings derived 
from the maps of ISchlegel et all d!998l) . 



5. Photometric zero points and absolute flux 
calibrations for use with the IRFM 

Eq. |4] is what needs to be implemented practically to derive 
IRFM temperatures. The quantities on the lefthand side are ob- 
served quantities while those on the righthand side are theoreti- 
cal quantities. One is immediately faced with a series of choices 

1 . the magnitude of the standard star (m ca i) 

2. the monochromatic flux of the standard star (F c . d \) 

3. the zero point for q 

4. the zero point for R,i wo 



These choices are only apparently trivial. The 2MASS mag- 
nitudes have been carefully calibrated in absolute fluxes by 
ICohen et alJ (|2003); however, the standard star to which the 
whole system is tied, Vega, has not been observed by 2MASS 
with sufficient accuracy due to its high brightness. A possible 
solution is to assume that the 2MASS magnitudes of Vega are 
given by the zero points of Cohen et al., with changed sign, as in 
ICasagrande et a l. (2006). Another complication is the zero point 
of the theoretical quantities. It is obvious from the definition of q 
that its value is 1 for the standard star; however, what are the cor- 
rect effective temperature, metallicity, and surface gravity of the 
standard star? The Cohen et al. calibration relies on an ATLAS 
theoretical spectrum computed by R.L. Kurucz with the "OLD" 
opacity distribution functions assuming T e s = 9400, \ogg = 3.9, 
a metallicity of -0.5, and a microturbulent velocity of km s _1 . 
Such a spectrum is not available in tabular form, we could indeed 
recompute it, however using such a spectrum to zero our theo- 
retical quantities would mean using a spectrum that is computed 
from a model inconsistent with the rest of the theoretical grid. 
Furthermore, as we shall see in Sect. [6] we will also need the ab- 
solute fluxes in other bands to derive a calibration for the bolo- 
metric flux. The natur al choice is to use th e corresponding the- 
oretical magnitudes o f lBessell et al.1 (jl998) transformed into the 
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Fig. 3. The same as Fig. [2] but for giants. 



Table 1. Monochromatic Fluxes for Vega from the calibrated 
ATLAS 9 flux. 



Band 


Wavelength 


Flux 


Mag Vega 




nm 


10-9 








erg s -1 cirT 2 nm"' 




J 


1235 


3.072 


0.038 


H 


1662 


1.113 


0.040 


K 


2159 


0.418 


0.043 



2MASS syst em. These magnitudes rel y on the model for Vega, 
pr oposed b y Castelli & Kurucz ( 1994), consistent with the grid 
of ICastelli & Kurucz I (12003^ that we are using. A possible solu- 
tion is to follow what was done by Casagran de et all (12006). who 
in fact used two different calibrations for optical and IR mag- 
nitudes. Inspection of Eq. [4] suggests another solution: use the 
same spectrum of the standard star to calibrate all bands. In this 
way any error in the calibration will cancel out when computing 
the flux ratio on the lefthand side of Eq. [4] However, to have a 
good absolute calibration, one also needs accurate observed or 
derived 2MASS magnitudes of the standard star Vega, which is 
quite difficult to obtain. We decid ed to adopt as 2 MASS mag- 
nitudes of Vega those provided by iMcCall ( 120041 see Sect. 
This theoretical spectrum should also be used to define the zero 
point of q and /? t heo for the standard star, in order to have a fully 
self-consistent IRFM. 



Throughout this work w e adopt the theoretical flux of Vega 
of [Castelli & K uruczl (fT994)l 5 l which has been calibrated to ab- 
solute flux, at Earth, using the value recommended by lHavesI 
(I1985L 3.44 x 10~ 8 erg s _1 cirT 2 nm -1 ). This spectrum is used 
to define the zero point of the q factor and the monochromatic 
fluxes, at the isophotal wavelengths of the 2MASS filters, listed 
in Table [H are used in ou r implementation of Eq. [4] As noted by 
Casagran de et al.l d2006l) . such a calibrated spectru m differs to 
the on e used by Cohen et al. (1992) and adopted by Cohe n et al.l 
d2003l) to define the absolute flux calibration of the 2MASS 
magnitudes. In Fig. [T] we display the calibrated spectrum of 
Vega in comparison with the adopted monochromatic fluxes of 
Cohen et al. (2003). The difference is small when comparing 
them with the continuum flux of our ATLAS 9 model of Vega 
at the same infrared wavelenghts. We stress that, for the purpose 
of consistent IRFM temperatures, we are not all that interested 
in having accurate monochromatic fluxes, but instead accurate 
ratios of bolometric fluxes to monochromatic fluxes. However, it 
should be noted that the adopted observed magnitudes for Vega 
are perhaps the source of uncertainty in a given temperature scale 
based on the IRFM. 

We adopted an error of 1 per cent on the monochromatic flux 
of each band for the determination of effective temperatures. The 
average systematic errors due to this uncertainty are 44 K and 
46 K for dwarfs and giants, respectively. These errors were es- 



5 http://wwwuser.oats.inaf.it/castelli/vega/fm05t9550g395k2odfnew.dat 
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Table 3. Adopted absolute integrated fluxes and magnitudes for 
Vega. 



Band 


Flux 


mag(Vega) 




10"5 






-1 -2 

erg s cm 




U 


0.267 


0.024 


B 


0.607 


0.028 


V 


0.321 


0.030 


R 


0.341 


0.037 


1 


0.167 


0.033 


J 


0.050 


0.038 


H 


0.028 


0.040 


K 


0.011 


0.043 



timated by adding the errors quadratically on temperature from 
each band and calculating the average over all stars in both sam- 
ples. 

6. Bolometric fluxes 

One of the fundamental observational quantities for applying the 
IRFM is the bolometric flux. The bolometric flux is not read- 
ily available for any given star, but Bl ackwell & Petfordl dl99ll) 
suggested that one could use a relation of the type F\, \ = 
10~°- 0(X, [Fe/H]), where m is a suitable broad band colour 
and X a colour ind ex. Such a calibration has been derived by 
lAlonso et all d 1995b using the K magnitude and the V — K colour 



Alonso et al. ( using the K magnitude ana tne V - K colour 
(Johnson system); Casagrand e et alj ( 120061) have derived several 
similar calibrations for different choices of m and the colour in- 
dex. In an initial attempt we tried to use the lAlonso etaild 1995b 
calibration for this purpose, which provided satisfactory results; 
however, the referees have correctly pointed out that, in doing so, 
we were forced to transform our (V - K s ) colour into Johnson's 
system, thus losing the internal consis tency of the 2MASS sys- 
tem. Furthermore, we had to apply the Alo nso et al.l (Tl995) cali- 
bration outside its formal ran ge of applicability, f or very metal- 
poor stars. The calibration of 1 Alonso et al.1 (1 19951) only had two 
stars at [Fe/H] = -3.2 and -2.9 and the rest with [Fe/H] > -2.6. 
A similar extrapolation proble m would apply if w e had u sed any 
of the calibrations derived by Casagr ande et al.l d2006l) . which 
were derived for stars with [Fe/H] > -1.9. We therefore decided 
to derive a new calibration that makes use of the K s magnitude 
and the 2MASS-based (V-K s ) colour and covers the metallicity 
range appropriate to our sample of stars. 

We adopt an approach si milar to that of lAlonso et al.l d 19951) 
andlCasagrand eetal] d2006). with a slight difference. The above 
authors use a set of effective wavelengths and monochromatic 
fluxes for Vega in order to define the integrated flux within each 
broad band from the photometry and the magnitudes of Vega. 
From the definition of magnitude follows 

F„ = F VegR 10- QA( - m - m ^ (5) 

where m is any photometric band. Provided then that the inte- 
grated flux of Vega in any given band is known, the integrated 
flux for the target star may be simply derived from its measured 
magnitude and the magnitude of Vega. In Table [3] we provide 
our adopted integrated magnitudes for Vega for the bands we 
are interested in, and JHK S refer to the 2MASS colours. These 
integrated magnitud es were d erive d by integrating the filter re- 
sponse func t ions o f iBesselfl d 19901) for the optical bands and 
ICohen et all d2003l) for the 2MASS bands, over the theoreti- 
cal flux of Vega. Consistently, the magnitudes for Vega were 



ta ken from McCall (2004), which gives the model magnitudes 
of iBessell et al.l d!998l) for the optical bands. This author cal- 
culates the 2MASS magnitu des of Vega using t he IR absolute 
monochromatic fluxes from Coh en et ail d2003l) . which are in 
fact quite similar to our adopted absolute monochromatic fluxes 
from the calibrated model of Vega (see Fig. [T). Initially we 
were going to adopt the magnitudes of Vega equal to zero in 
the 2MASS bands, but when we derived the IRFM tempera- 
tures, our scale of te mperatures was ~ 120K hotter than that of 
lAlonso et al.l d!996al) . which we think is the best implementation 
of the IRFM available in the literature due to its internal consis- 
tency in the whole range of metallicities f rom -3.0 to 0. 5. By 
adopting the magnitudes of Vega given bv lMcCalll dlOOl . this 
difference is reduced to ~ 60 K, which we consider more appro- 
priate (see Sect. I8.lt . In addition, the remarks on zero points of 
Sect.|5]also apply here. This choice guarantees that any error in 
the absolute calibration of bolometric fluxes and monchromatic 
fluxes will cancel out in Eq. [4] 

From Eq. and the data in Table [3] for any star for which 
photometry in several bands is available, one may computed the 
total flux measured at Earth in the given bands. From this value, 
the bolometric flux can be obtained by using model data. 

In practice the stars we decided to use to derive the bolomet- 
ric flux calibrations fall into one of four groups: 

1. stars with full UBV(RI) C JHK S data, 

2. stars with UBVJHK S data, 

3. stars with BVJHK S data, 

4. stars with BV(RI) C JHK S data. 

The information is more complete for the stars of group 1) 
and 4); however, we decided to include the also the stars of 
groups 2) and 3) in our calibration effort, since this allows us to 
include a relevant number of stars at extremely low metallicity. 

For each group one may compute 



obs 



F(X)Ti(A)dA 



(6) 



where the sum is extended to all the available bands and T,(A) 
is the response function of the i - th band, the integrals must 
be considered extended from to infinity, formally, although in 
practise the response function of each filter vanishes outside a 
finite interval and numerically one stops integrating outside this 
interval. One can then use the models to compute the correction 



F boi - AF b s 



(7) 



where A is a function of T e g , logg and [Fe/H]. Obviously 
a different A has to be computed for any given set of available 
bands. The A factors for the different band combinations we have 
used are given in machine readable form at the CDS (see Sect.O 
with different subscripts 1 to 4 corresponding to the different 
band combina tions. This i s again sli ghtly different from wha t 
was done by lAlonso et ail d 19951) or ICasagrande et alj d2006), 
who used the monochromatic fluxes at the effective wavelengths 
of each band to approximate the spectral energy distribution of 
the star and integrated this approximate energy distribution over 
the whole interval. We only make use of integrated fluxes, which 
are related to the observed magnitudes through Eq.|5]and of the 
fluxes and magnitudes of Vega given in Table[3] 

While the 2MASS magnitudes are provided by the cata- 
logue, it is customary for the optical bands to provide the V 
magnitude and the colours U - B, B — V, V — R, V — I. While 
for the bands B to / it is straightforward to obtain the magnitude 
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Colour 


Colour range 


[Fe/H] range 


a" 
"o 


a\ «2 


"3 


04 


"5 


«6 




o>(%) c 


Dwarf stars 


V- 


J 


[0.8,2.4] 


[-3.5,0.3] 


2.4945 


-2.2635 0.9615 


-0.1509 


0.0657 


-0.1365 


-0.0074 


219 


0.9 


V - 


H 


[0.9,3-0] 


[-3.5,0.3] 


2.3681 


-1.9055 0.6415 


-0.0773 


0.0418 


-0.1028 


-0.0053 


216 


0.9 


V - 


a; 


[1.0,3.0] 


[-3.5,0.3] 


2.3522 


-1.8817 0.6229 


-0.0745 


0.0371 


-0.0990 


-0.0052 


216 


0.9 


Giant stars 


V- 


j 


[0.5,2.7] 


[-4.0,0.1] 


2.2282 


-1.7818 0.6809 


-0.0923 


0.0302 


-0.0696 


-0.0031 


97 


1.3 


V- 


H 


[0.6,3.4] 


[-4.0,0.1] 


2.1522 


-1.5792 0.4821 


-0.0523 


0.0182 


0.0502 


-0.0019 


91 


1.3 


V - 


a; 


[0.7,3.8] 


[-4.0,0.1] 


2.1304 


-1.5438 0.4562 


-0.0483 


0.0132 


0.0456 


-0.0026 


95 


1.4 



" The coefficients of the calibrations a, are given in units of 10~ 5 erg cnT 2 s _I . 

* N is the remaining number of stars after several iterations (usually less than 20) of the 2.5cr clipping. 
cr 0, given in per cent, is the standard deviation of the final calibrations. 
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Fig. 4. Comparison of our temperature scale with that of lAlonso et al.1 {1996a, 1999a). The dashed-dotted line indicates the average 
temperature difference and dashed lines the standard deviation, l<x, from the average (see text). 



(e.g. R - V - (V - R)), some caution must be exerted for the 
U band, for which the atmospheric extinction is strongly vary- 
ing across the band and establishes the UV cut-off. In f act, this 
band has proved to be the most difficult to standardize. iBesselll 
(1990) provides a response curve UX in which the atmospheric 
extinction is folded in and gives the curves BX and B for the B 
band. The former is to be used to compute the synthetic pho- 



tometry of the (U - B) colour, while the latter is to be used to 
compute the synthetic (B - V) colour. For the stars to be used in 
our bolometric flux calibration we define the U magnitude 

U = (U-B) + V + (B-V) + (BX-B) (8) 

where V, (U - B), and (B - V) are the observed magnitude and 
colours of the star, and (BX - B) is derived from the theoretical 
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Fig. 5. Left: Empirical calibration T e ff:(B - V)-[Fe/H] for dwarfs in the metallicity bins -0.5 < [Fe/H] < 0.5 (filled circles), 
—1.5 < [Fe/H] < -0.5 {open circles), -2.5 < [Fe/H] < —1.5 (triangles), and [Fe/H] < -2.5 (diamonds). The lines correspond to 
our calibration for [Fe/H] = (solid line), -1.0 (dotted line), -2.0 (dashed line), -3.0 (dotted-dashed line). Right: Residuals of the 
fit (Ar eff = T^ FM - T C J) as a function of (B - V). 



models with the requirement that it must be equal to zero for 
Vega. The colours (BX - B) are given in the online data at the 
CDS (see Sect. |3). 

The bolometric flux of each star was determined as in 
lAlonso et all 0995). We first determined the fluxes of each band 
by applying Eq.|5]from the observed magnitudes of the star. Then 
we derived the bolometric flux, Fb i, using the Eq. |7j Thus the 
temperature T e g was then determined using the IRFM that com- 
bines fboi and the monochromatic fluxes at IR wavelengths. This 
new value for the effective temperature may re-enter in Eq.UJto 
derive a new value for F\, \, and so on. This iterative procedure 
converges quickly towards a final Fboi- In each iteration, the fac- 
tors A and BX-B were determined using a trilinear interpolation 
within the grid for the corresponding T e ff , log g, and [Fe/H] of 
the star. We considered the errors on the IRFM T e ff due to uncer- 
tainties on the adopted absolute calibration of the 2MASS pho- 
tometric system, and the errors on the magnitudes JHK^ log g, 
and [Fe/H]. 

We derived relations between bolometric fluxes and colours, 
taking also the effects of metallicity i nto account. We adopte d 
the same fitting formula as adopted by Casagrande et al. (2006) 



[Fe/H]) = a + a x X + a 2 X 2 + arf 1 + 

+ a 4 X[Fe/H] + a 5 [Fe/H] + a 6 [Fe/H] 2 (9) 

where the cf> is derived as <p(V - m, [Fe/H]) = F bo i/10-°- 4m , X = 
V - m represents the J,H,K S magnitudes, and a,- (z = 0, ...,6) 
are the coefficients of the fit. We iterate the fitting procedure by 
discarding the points more than 2.5<r from t he mean fit. We als o 
tried other fitting formulae such as that of lAlonso et al.1 ( 1 1995b . 
but they led to similar re sults. This has been extensively tested by 
ICasagrande et al.1 (12006). who also give fits using optical bands. 

For dwarf s tars, we selected all the s tars of group 1) from 
the sam ple of lAlonso et al.1 dl996al) and [Ramirez & Mel endez] 
(2005a) with uncertainties in the JHK S magnitudes less than 0.1 
at all metallicities. We added t he stars of group 4) from the sam- 
ple o flCasagrande et aLl (120061) with the same restrictions and we 
completed the sample with stars of groups 2) and 3) with [Fe/H] 
< -2. We only added metal-poor stars to give more weight to the 
metal-rich stars of groups 1) and 4); otherwise, we would end up 
with including all dwarf stars and the fit would be dominated by 
the greater number of stars of group 3). 

For giant sta rs, we applied the sam e re strictions but the num- 
ber of stars with lAlonso et al.1 (Il996al) and lRarmrez & Melendez! 
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Fig. 6. The same as Fig. [5] but for V - K s in dwarfs. 



(2005i of groups 1) and 4) was very small (only 10 stars in 
group 1) and 26 in group 4) with uncertainties in the JHK S mag- 
nitudes less than 0.1), so we decided to include all the stars of 
groups 2) and 3). 

In Figs. [2] and [3] we display the polynomial fits that repre- 
sent the empirical calibrations (f> versus colours and metallicity. 
The coefficients of these calibrations are given in Table [2] to- 
gether with the remaining number of stars after the 2.5<x clip- 
ping and the r.m.s. of the fit, cr^. These calibr atio ns show simi- 
lar be haviours to those of lAlonso et all (1 19951) and lAlonso et al.1 
(1999a), at least in the metallicity range from -3 to 0. 



7. IRFM temperatures and angular diameters 

To determine effective temperatures we need to apply Eq.|4] The 
bolometric fluxes are estimated using the empirical calibration 
F ho \M = lQ-° AK *<t>(V - K„[Fe/H]) given in Tableland the 
2MASS K s and Johnson V magnitudes. The q- and /^-factors 
are determined from an initial guess of the temperature of the 
star, r e fFo, by trilinear interpolation in the grid, using the sur- 
face gravity and metallicity of the star. Then, we determine a 
new value for the temperature by comparing the theoretical bolo- 
metric flux, FboUheo, derived from the previous determination 
of T e ff and the bolometric flux, Fboi.cah using the expression: 
7 , efr,new =T ef f , id[^boi,cai/^ Boi.theo] 1/4 ■ We again derive the q- and 



/^-factors for r e (f new and repeat this process iteratively until 

|7eff,new-r e ff,oldl < 0.1 K. 

The final temperature of the star is determined as the aver- 
age of the three temperatures extracted from each of 2MASS 
filters weighted with the inverse of their individual errors (see 
lAlonso et al.lll996ab . The error on the weighted mean is com- 
puted as Ar e ff - Nl Yj (Ar,-) -1 where Ar,- are the errors of the 
temperatures from the individual filters (i = J,H,K S ) and N = 3 
is the number of available temperatures. These errors Ar, ac- 
count for the photometric errors of the observed JHK S and V 
magnitudes, the error on the adopted absolute calibration for 
the 2MASS photometric system, and the uncertainties on sur- 
face gravity and metallicity. To estimate Ar,-, we just add all the 
individual errors of the i band quadratically. 

The angular diameters have been calculated from Eq.[TJwith 
the derived IRFM temperatures and bolometric fluxes. Their er- 
rors were estimated by propagating a mean error of 1.3% in the 
bolometric fluxes and the errors on T e ff . 



8. Comparison with other temperature scales. 

In this section, we compare our temperature scale with other 
temperature determinations based on different implementations 
of the IRFM (lAlonso et al.lll9961j[r999atlRarmrez & Melendezl 
12005 a; Casagrande et al. 2006), on the excitation equilibrium of 
Fei lines (San tos et alj |2004). and on the fitting of Balmer line 
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Table 4. Comparison with other temperature scales. 



Sample 


[Fe/H] range Ar cff a TdY 


N" 


Dwarf stars 



Alonso et al. (1996a) 
Alonso et al. (1996a) 
Alonso et al. (1996a) 
Ramirez & Melendez (2005a) 
Ramirez & Melendez (2005 a) 
Ramirez & Melendez (2005 a) 
Casagrande et ah (2006) 
Casagrande et al. ( 2006)* 
Santos et al. (2004) 
Sant os et al. (2004) * 
Boni facio et al. (2007) 
Barkl em et al. (2002) 
Bark lem et al. (2002) 
Christlieb et al. (2004) 
Baines et al. (2008) c 



r_s 5 +0 SI 

L J.J) T U.JJ 


+64 


104 


332 


[-3.5,-2.5] 


+61 


91 


18 


[_ Ut J^TU. JJ 


+32 


130 


122 


[-4A+0.3] 


+33 


98 


84 


[^t.0,-2.5] 


-87 


194 


12 


[-0.5,+0.3] 


+45 


91 


69 


[-1.9,+0.4] 


-12 


56 


101 


[-1.9,+0.4] 


-41 


50 


101 


[-0.7,+0.5] 


+ 11 


120 


133 


[-0.7,+0.5] 


-13 


129 


133 


[-3.6,-2.4] 


+165 


79 


16 


[-2.5,+0.1] 


+77 


133 


23 


[-0.5,-0.1] 


+51 


129 


16 


[-3.1,-1.6] 


+ 177 


80 


8 


[-0.4,0.5] 


-32 


163 


22 



Giant stars 



Alonso etal. (1999a) 
Alonso etal. (1999a) 
Alonso et al. (1999a) 
Ramirez & Melendez (2005a) 
Ramirez & Melendez (2005 a) 
Ramirez & Melendez (2005a) 
Cavrel et al. (2004) 
Christlieb et al. (2004) 
Baines et al. (2008) c 



of 



profiles dBarklem et al.l l2002; Bon ifacio et aT1l2007l) . In Table 
we gather the mean differences between our temperatures and 
those of different samples, Ar e ff , together with the standard de- 
viation (scatter) around the mean, <x rEff . 

8.1. Alonso etal. sample 

Th e updated temperatures do no t differ significantly from those 
of lAlonso et ail dl996a. 1999 j see Fig. B). O ur" temperature 
scale is hotter than that of lAlonso et all dl996al Il999al) for both 
dwarfs and giants. We find an average difference Ar e ff = +64 K 
with a cr TeB = 104 K (N = 332 dwarfs) and A7/ eff = +54 K 
with a cr TcS = 131 K (A^ = 202 giants). This translates into 
a mean T e g difference of < 1%. Although not negligible, 
such differences are within the error bars of the current tem- 
perature determinations, although the scatter, ctt m , seems to 
be quite large. The different bolometric flux calibration, pho- 
tometric data and absolute flux calibration might be responsi- 
ble for this small differe nce between the two temperature scales. 
Casagrande et al. (2006) checked that using the absolute calibra- 
tion of I Alonso et al.ld 1995b and if using the TCS fil ters, their cal- 
ibration and that of iRamfrez & Melendez! d2005al) agree within 
20 K. However, to do this exercise they had to transform the 



2MASS magnitudes into the TCS system, so their conclusions 
may be affected by these transformations. 

Even if we select subsamples of different mean metallicity, 
the differences remain very small (see Tabled. In conclusion, in 
the whole metallicity range, the systematic d i fference betwe en 
our temperature scale and that of Alon so et al J d 1 996al 1 1 999ah in 
dwarfs and giants is positive but smaller than +65 K, which is 
in fact less that the average of the individual uncertainties in our 
calibration (< 82 K for dwarfs and < 76 K for giants). 

8.2. Ramirez & Melendez sample 



IRamfrez & Melendez! d2005al) 
of lAlonso et all d 1 996al 1 1 999al) 
from 



[-3.0.+0.5] 


+54 


131 


202 


[-3.0,-2.5] 


+76 


120 


10 


[-0.5,+0.5] 


+43 


144 


116 


[-4.0,+0.3] 


+63 


57 


25 


[-4.0,-2.5] 


+61 


62 


18 


+0.2 d 


+ 116 




1 


[-4.0,-2.0] 


+ 115 


76 


34 


[-3.4,-2.6] 


+ 128 


71 


22 


e 


-67 


139 


6 



The number of stars. 

If we consider all reddening corrections equal to zero. 
Ar c ff = r'^ FM - rij, where is a direct determination 
r eff using the angular diameter 6. 
One metal-rich giant star. 

Did not find any metallicity determination so decided to adopt 
[Fe/H] = 0. 



extend the sample of stars 
with metal-rich stars 
Sant os et al .1 d2004l) an d very meta l- poor stars from 
IChristlieb et all (120041) and ICavrel et all ((20041) . We de- 
termined effective t empera tures for the calibrators of 
IRamfrez & Melendez! d2005al) using our implementation of 
the IR FM. Our temperature scale i s on average hotter than 
that of IRamfrez & Melendez] d2005al) by about A7/ eff = +33 K 
(o-j eff = 98 K, N = 84 dwarfs) and Ar eff = +63 K (cr re(r = 57 K, 
N - 25 giants). This difference might be partially related to the 
use of different absolute calibration as we stated in Sect. 18.11 

Among giants, we find minor differences when we look at 
the most metal-poor and metal-rich stars in the sample (see 
Table |4j. However, this behaviour changes when we inspect the 
dwarf stars. While for metal-rich dwarfs we find Ar e ff = +45 K 
(o- TcB = 76 K, N = 69 dwarfs with [Fe/H] > -0.5), for metal- 
poor dwarfs we find our temperature scale to be cooler: AT e ff = 
-87 K (o- Tea = 194 K, N = 12 dwarfs with [Fe/H] < -2.5). 
We believe that this difference is mainly due to the photometric 
transformations between the 2M ASS and the TCS systems that 
Ramirez & Melendez (2005a) need to perform in order to derive 
the IRFM temperatures. 

8.3. Casagrande et al. sample 

Casagran de et al.l (|2()06) propose a new IRFM using multiband 
photometry. They derive empirical effective temperature and 
bolometric flux calibration for G and K dwarfs stars in the 
range -1.87 <[Fe/H] < 0.34. They use BV(RI) C Johnson- 
Cousins photometry and JHK S 2MASS photometry. We ap- 
plied our implementation to estimating the r e ff of stars in their 
sample and find our temperature scale only slightly cooler by 
Ar eff = -12K (cr Tca = 56 K, N = 101 dwarfs). For these stars, 
w e estimated the redde ning corrections from the maps of dust 
of ISchlegel et all dl998l) . corrected as described in Sect. 03] In 
Table [4] we also comp are our temperature scale with that of 
Casagrande et al. (2006) by arbitrarily adopting the reddening 
corrections E(B - V) = for all the stars, under the assumption 
that all these stars are nearby and should not show any reddening 
at all. We find AT eff = -41 K (o- TeS = 50 K, N = 101 dwarfs), 
i.e. temperatures 29 K cooler on average. This systematic differ- 
ence is well within our error budget, so we decided to use these 
reddening corrections to be consistent with other new stars in- 
cluded in the sample for which we need to estimate the redden- 
ing corrections as the dwarf stars from Bonifac io et al.l d2007f) 
(see Sect.lOll. 

8.4. Santos et al. sample 

ISantos et al.l d2004l) have carried out a detailed spectroscopic 
analysis of planet-host stars and a comparison sample of stars 
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Fig. 7. The same as Fig. [5] but for giants. 



without known planets. Their effective temperatures are based 
on the excitation equilibrium of the F e i lines. Our t emper ature 
scale is only slightly hotter than that of Sant os etldl {2004) with 
Ar eff = +11 K (tr^ = 120 K, N = 133 dwarfs), although with 
a large scatter. As in the previ ous section, we also derived the 
reddening corrections from the Schleg eTet al.l JT998) dust maps, 
corrected as described in Section 1431 In Table |4] we also show 
the comparison with E(B - V) = 0. 



8.5. Cayrel etal. sample 



8.6. Bonifacio et al. sample 



ICavrel et all (120041) present UVES spectroscopic observations 
of very metal-poor giant stars. T hey derived T ^g usin g the 
r eff :colour-[Fe/H] calibrations of lAlonso et all (Il999bh . Our 
r etf scale is hotter by AT eff = +115K (o- TeB = 76 K, N = 34 
giants with [Fe/H] < -2.5); however, that part of this difference 
stems from the different choice made for the re ddening. Here 
we ad opted the reddening from the dust maps of Schlegel et al. 
( 1998), corrected as de scribed in Sect.|4.3l where as ICavrel et all 
(2004) instead used the lBurstein & Heilesldl982l) maps. The dif- 
ferent choice in reddening account s for a difference of ~ 40 - 
50 K, on average (ICavrel et al.l2004l) . The remainin g 75 K reflect 
the diff erence between our calibration and that of lAlonso et al.1 
(1999b). Unsurprisingly, this is, essentially, the same as what 
was found for giant stars with [Fe/H] < -2.5 in Sect. 18. ll (see 
Table©. 



Boni facio et alJ(l2007l) present high quality spectroscopic data of 
a sample of extremely metal-poor dwarf stars. They derived the 
effective temperatures by fitting the wings of the Ha line. We de- 
rived the effective temperatures of these stars using the 2 MAS S 
JHK S magnitudes and reddenings from Schleg el et al.l l[l998), 
corrected as described in Sect. 14.31 Our effective temperatures 
are significantly hotter than those derived from the Balmer lines, 
Ar eff = +165K (ct Tm = 79 K, N = 16 dwarfs with [Fe/H] < 
-2.5). The difference between the temperatures derived from 
Ha, an d those derived by using the co lour V - K in the calibra- 
tions of IRamfrez & Melendez)l2005bl) is roughly 265 + 122 K. 
This ~ 100K difference may be partiall y explained by our 
compa rison with the temperature scale of IRamfrez & Melendezl 
(l2005al) with Ar eff = -87 K (see Sect. This has an im- 
pact on the Li abun dances in extremely m etal-poor stars down 
to [Fe/H] = -4 (see Sbo rdone et afll2008l in prep.), because the 
IRFM temperatures would provide higher Li abundances at the 
lowest metallicities, whereas Ha temperatures seem to show a 
slowly decreasing trend in Li towards lower metallicities. 

iBo nifacio et al. (2007]) use the theory of Ba rklem et al.l 
(2000) to describe the self-broadening o f Balm er lines. For the 
same sample of stars, IBonifacio et al.l d2003l) instead use the 
lAli & Grieml d!965l 1 19661) theory and derived effective temper- 
atures which were on average 150 K hotter, thus in substantial 
agreement with our IRFM temperatures. 
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8. 7. Christlieb et al. sample 

IChristlieb et all (|2004|) present the Hamburg/ESO R-process 
Enhanced Star survey (HERES) with the aim of searching for 
very metal-poor stars ([Fe/H] < -2.5) with r-process elements 
enhanced. We selected those s t ars wi th available B — V and V 
photometry in Christlie b et al.l (120041) and took the stellar pa- 
rameters from Bark lem et al.l d2005l) . The effective temperatures 
were estimated by averaging the resultin g T e g from the dif - 
fere nt r e ff:colour - [Fe/H] calibrations of lAlonso et alj (Q996b) 
and lAlonso et al.l dl999bl). They followed the prescription de- 
scribed by ISivarani et al.l (120041) . Our T e g scale is significantly 
hotter by AT eff = +177K (cr rdT = 80 K, N = 8 dwarfs with 
-3.1 [Fe/H] < -1.6) and Ar eff = +128 K (o- Tclr = 71 K, N = 22 
giants with -3.4 [Fe/H] < -2.6). This difference probably comes 
from the different adopted reddenings a nd the difference be- 
twe en our temperature s cale and that of lAlonso et al.l ([T996a) 
and lAlonso et all dl999al) . 



effects add to the systematic error budget of Balmer-line based 
effectiv e temperatures . We co mpared our temperature scale with 
that of Barkle nTet alj (120021) . Our temperatures are hotter by 
Ar eff = +77 K {cr TtS = 133 K, N = 23 dwarfs with [Fe/H] > 
-2.5). However, for metal-rich dwarfs with [Fe/H] > -0.5, this 
difference drops to AT eff = +5 1 K (cr TdT = 129 K, N = 16 dwarfs 
with [Fe/H] > -0.5). The average difference between the two 
temperature scales remains within the uncertainties on the tem- 
perature determinations, although the standard deviation is large. 



9. r eff :colour-[Fe/H] calibrations 

We derived relations between T e ff and colours, also taking the 
effects of metallicity into accou nt. We adopt e d the same fit - 
ting formula that was adopted bv lAlonso et all dl996blfl 999b). 
iRarrurez & Melend ez (2005bj), and lCasagrande et all d2006l) 



8.8. Barklem etal. sample 

Balmer-line profile fitting in principle allows a very precise 
determination of ste llar effective temperature for cool stars. 
Bark lem et al.l d2002l) claim an accuracy of the temperature de- 
terminations of ~ 65 K for solar metallicity stars but for [Fe/H] ~ 
-1 of ~ 80 and [Fe/H] ~ -2 of ~ 100 K. Uncertainties in the 
theory of self-broadening, deviations from LTE and granulation 



eff = b +biX+b2X 2 +b 3 X[Fe/U]+b4[Fe/H]+b 5 [Fe/H] 2 (10) 

where 6 e ff = 5040/ T e s, X represents the colour, and b, (i = 
0, 5) are the coefficients of the fit. We iterate the fitting proce- 
dure by discarding the points more than 2.5<x from the mean 
fit. All our calibrations were adequately tested by removing 
some terms and/or adding higher order terms in either X and 
[Fe/H] . We verified that neither removing terms nor introducing 
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Table 5. Coefficients and range of applicability of the T e ff :colour-[Fe/H] calibrations. 
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Colour 


Colour range 


[Fe/H] range 


bo 


bi 


b 2 


b 3 


b 4 


b 5 


N" 


°X ff (K) 6 


Dwarf stars 


B- 


V 


[0.2,1.3] 


[-3.5,0.5] 


0.5725 


0.4722 


0.0086 


-0.0628 


-0.0038 


-0.0051 


418 


76 


V- 


R 


[0.2,0.8] 


[-3.1,0.3] 


0.4451 


1.4561 


-0.6893 


-0.0944 


0.0161 


-0.0038 


164 


45 


V- 


I 


[0.5,1.4] 


[-3.1,0.3] 


0.4025 


0.8324 


-0.2041 


-0.0555 


0.0410 


-0.0003 


164 


52 


V- 


J 


[0.5,2.3] 


[-3.5,0.5] 


0.4997 


0.3504 


-0.0230 


-0.0295 


0.0468 


0.0037 


430 


36 


V- 


H 


[0.6,2.8] 


[-3.5,0.5] 


0.5341 


0.2517 


-0.0100 


-0.0236 


0.0523 


0.0044 


426 


30 


V- 


K s 


[0.7,3.0] 


[-3.5,0.5] 


0.5201 


0.2511 


-0.0118 


-0.0186 


0.0408 


0.0033 


431 


32 


J - 


K s 


[0.1,0.8] 


[-3.5,0.5] 


0.6524 


0.5813 


0.1225 


-0.0646 


0.0370 


0.0016 


436 


139 


Giant stars 


B- 


V 


[0.3,1.4] 


[-4.0,0.2] 


0.4967 


0.7260 


-0.1563 


0.0255 


-0.0585 


-0.0061 


120 


57 


V- 


R 


[0.3,0.7] 


[-4.0,0.1] 


0.4530 


1.4347 


-0.5883 


-0.0156 


-0.0096 


-0.0039 


55 


85 


V- 


J 


[1.0,2.4] 


[-4.0,0.2] 


0.4629 


0.4124 


-0.0417 


-0.0012 


0.0094 


0.0013 


138 


18 


V- 


H 


[0.8,3.1] 


[-4.0,0.2] 


0.5321 


0.2649 


-0.0146 


-0.0069 


0.0211 


0.0009 


144 


23 


V- 


K s 


[1.1,3.4] 


[-4.0,0.2] 


0.5293 


0.2489 


-0.0119 


-0.0042 


0.0135 


0.0010 


145 


23 


J- 


K s 


[0.1,0.9] 


[-4.0,0.2] 


0.6517 


0.6312 


0.0168 


-0.0381 


0.0256 


0.0013 


145 


94 



" The remaining number of stars after several iterations (usually less than 20) of the 2.5<r clipping. 
* o"r cir , given in K, is the standard deviation of the final calibrations. 



higher order terms improves the accuracy of the fit significantly. 
Therefore we adopted Eq. [TO] 

In Figs. [5] |6]|2l and[8] we display the polynomial fits which 
represent the empirical calibrations T e $ versus colours and 
metallicity. We discarded all the stars with uncertainties in the 
JHKs magnitudes greater than 0. 1 . The coefficient of these cal- 
ibrations are given in Table [5] along with the remaining num- 
ber of stars after the 2.5cr clipping and the r.m.s. of the fit, 
0"(7eff). Normally, the number of iterations were fewer than 20. 
Our polynomial fits of the colour B - V usually have si milar 
r.m.s. than those provided by IRamfrez & Melendez! (l2005bl) . We 
should point out the small number of giant stars with metallici- 
ti es [Fe/H] > - 1 .5, bec ause most of the giant stars in the sample 
of lAlonsoetai] (fl999a) are very bright objects, hence with poor- 
quality 2MASS JHK S magnitudes. For dwarf stars, our empiri- 
ca l calibrations of V - R and V - 1 have smaller r.m.s. than those 
of Ramirez & Melendez (2005b). For gia nt stars, the calibration 
of V - R shows a greater r.m.s. than in [Ramirez & Mel endez] 
(2005b), probably due to the small number of stars in our sam- 
ple. 

On the other hand, our empirical calibrations of the colours 
V - J,V-H, and V—K s have a sma ller r.m.s. than those presented 
by iRamfrez & Melendez! d2005bl) . For giants, our fits are more 
accurate although o ur sample contains sl i ghtly fe wer giant stars 
than the sample of Ramirez & Mel endez] (l2005bl) . 

10. Angular diameters 

The IRFM was developed to provide the T e ff and 9 simulta- 
neously from observed and theoretical data. One fundamental 
test to the IRFM is thus the comparison with measured angular 
diameters. IRamfrez & Melendez! (l2005al) compare their results 
with m easu red angular di a meters from iRichichi & Percheronl 
d2002l) and iKervella et ail d2004 for giant and dwarf stars, 
respectively. These stars are too bright for the 2MASS cat- 
alogue, making the 2MASS JHK S magnitudes very uncer- 
tain. However, iBaines et al.l d2008l) have recently presented 
new measurements of angular diameters using the Center for 
High Angular Resolution Astronomy (CHARA) Array, a six- 




0.2 0.4 0.6 0.8 -0.4 -0.2 0.0 0.2 0.4 0.6 

6 L „ (mas) [Fe/H] 



Fig. 9. Compariso n between the angular diameters from 
Bain es et"aT1 (120081) . #ld> and diameters from the IRFM, 6>irfm 
for dwarf stars. Error bars are individual uncertaintes. The dif- 
ferences between #irfm and 6*ld, A0, are also shown containing 
the uncertainties of both 6>irfm and Old- 



element Y-shaped interferometric array. We searched f or the 
JHK S magnitudes of the stars reported in Bai nes etldl (120081) 
and the results are presented i n Figs. [9] and JTOl The stellar 
parameters were adopted from Bai nes et al.l (120081) . and the 
metallicities were extracted fr om ISantos et al.l (|2004) and from 
Cavrel de Stro bel et al.l (12001 1) . However, for same cases, espe- 
cially for giant stars, we did not find any available metallicity 
determination, so we decided to adopt [Fe/H] = 0. The V mag - 
nitudes were extracted from the GCPD (Mermill iod et al.fl997l) . 
and in those cases where no value was found, we took the V mag- 
nitude as given in the SIMBAD catalogue. One can compare our 
IRFM angular diameters, 6>irfm, with the direct measurements, 
6>ld- For dwarf stars, the average difference, A6 - (^rfm - Old 
is 0.002 with a standard deviation erg = 0.033 (N = 22 stars). 
We can also derive a direct temperature, T^, from ^ld and 
the bolometric flux, determined from our bolometric flux cali- 
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Fig. 10. Same as Fig. [9] but for giants. 



bration, by using Eq. Q] The previous comparison between an- 
gular diameters thus translates into a temperature difference, 



Ar eff = T 



IRFM 

eff 



- T™, of -32 with o- TclT = 163 (N = 22 stars). For 
giants, the number of stars with relative ly accurate JHK S data is 
low. The sample of Baines et al. (2008) contains only six giant 
stars. For these stars, we find A9 = 0.012 with a cr g — 0.029 
which translates into Ar e ff = -67 with a ctt m = 139. These re- 
sults are also given in Table|4]in comparison with other tempera- 
ture determinations. Our new implementation of the IRFM pro- 
vides good results when comparing with direct measurements 
of angular diameters. The absence of any trend with metallicity 
in the residuals shown in Fig.|9]over almost 1 dex in metallicity 
is very encouraging. This suggests that the model atmospheres 
correctly model the variation of fluxes with metallicity. Since the 
metal-rich range is the most difficult for modelling the opacity, 
it is reasonable to expect that the models are also reliable at low 
metallicity. In other words, we do not expect that our tempera- 
ture scale has spurious trends with metallicity due to inadequate 
modelling of the stellar atmospheres. 



11. Summary 

We have made use of the IRFM to determine effective tempera- 
tures of 555 dwarf and subgiant field stars and of 264 giant field 
stars. Our implementation of the IRFM uses the 2MASS photo- 
metric system as a reference system to perform all the calcula- 
tions. We derived a bolometric flux calibration down to metal- 
licities of [Fe/H] = -3.5 for dwarfs and -4.0 for giants, as a 
function of the 2MASS magnitude, m, and the colour, V - m. 
We computed theoretical magnitudes by integrating the ATLAS 
models in the 2MASS JHK S filters. 

Our tempe rature scale is hotter than that of lAlonso et alj 
d!996al [19991) by ~ 64 K (<r T = 104 K, N = 332 dwarfs) and 
~ 54 K (it t = 131 K, N = 202 giants). Simil a r result s are found 
when comparing with Ramirez & Mel endezl d2005al) . However, 
interestingly, for dwarfs stars with [Fe/H] < -2.5, the compar- 
ison with the sample of iRamfrez & Melendezl d2005al) provides 
a difference of ~ -87 K (<r r = 1 94 K, N = 12 dwarfs) We be- 
lieve this difference is related to Ra mirez & Melendezl (|2005a) 
transforming the 2MASS JHK S magnitudes to the TCS photo- 
metric system to derive bolometric fluxes and effective tempera- 
tures for their calibrators at the lowest metallicities, whereas we 
determine the bolometric fluxes and effective temperatures in the 
2MASS photometric system. 



Our r e ff are hotter th an those estimated us ing Ha line pro- 
files by AT eS = +77 K dBarklem et alj 120021 a Tm = 133 K, 
N = 23 dwarfs with [Fe/H] > -2.5) and Ar eff = +165 K 
dBonifacio et alj I2007L cr Tclr = 79 K, N = 16 dwarfs with 
[Fe/H] < -2.5). This result has implications for the Li abun- 
dances for very metal-poor stars down to [Fe/H] = -4. Higher 
temperatures provide higher Li abundances; therefore, the drop 
of the Li abundances towards lower metallicities will cancel out, 
and the Li abundances would remain in a plateau. 

We derived T e s versus colour empiric al calibration s , which 
are co mp atible with those pres ent ed by [ Alonso et al] (fl 996b , 
19999), IRamfrez & Melendezl d2005bE ICasagrand e et alj 
(2006), within the quoted errors. For those who wish to use 
2MASS photometry to estimate effective temperatures for a 
wide range of metallicities, we recommend our calibration, 
which were derived within the 2MASS system, rather than 
the others, which are either based on different systems or on 
hybrid systems. Our calibrations exploit the excellent internal 
consistency of the 2MASS photometry and should provide 
accurate temperatures in a relative sense. In an absolute sense, 
our calibrations are of the same quality as the other calibrations. 

A comparison of IRFM angular diamete rs with interferomet- 
ric measurements of angular diameters from Bai nes et al. (2008) 
provides good agreement for both dwarf and giant stars. This 
gives us confidence that our new implementation of the IRFM is 
reasonable. 
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